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We propose an approach to quantum phase estimation that can attain precision near the Heisen¬ 
berg limit without requiring single-particle-resolved state detection. We show that the “one-axis 
twisting” interaction, well known for generating spin squeezing in atomic ensembles, can also amplify 
the output signal of an entanglement-enhanced interferometer to facilitate readout. Applying this 
interaction-based readout to oversqueezed, non-Gaussian states yields a Heisenberg scaling in phase 
sensitivity, which persists in the presence of detection noise as large as the quantum projection noise 
of an unentangled ensemble. Even in dissipative implementations—e.g., employing light-mediated 
interactions in an optical cavity or Rydberg dressing—the method significantly relaxes the detection 
resolution required for spectroscopy beyond the standard quantum limit. 


For decades, advances in atomic spectroscopy have 
brought clocks, accelerometers, and magnetometers to 
ever greater precision. A recent development is the use 
of many-particle entangled states to reduce the statistical 
uncertainty in measurements of the energy difference huo 
between two atomic states \f), \i) [THTO]. Whereas an un¬ 
correlated ensemble of N two-level atoms achieves at best 
the standard quantum limit of precision A uoT = 1/y/N 
in an interrogation time T, entanglement can enhance 
this precision up to the fundamental Heisenberg limit 
Ac oT = 1/N. Approaching the Heisenberg limit with 
more than a few particles remains a major outstanding 
challenge, due to difficulties not only of preparing but 
also of detecting entangled quantum states QTH33]. 

Imperfect state detection has limited the sensitivity of 
entanglement-enhanced metrology with squeezed w, 
oversqueezed j9], and twin Fock m states. The standard 
detection protocol is to measure the population difference 
n = ?7q — 77q between the levels |t), \i) in the entangled 
state after perturbing it by an amount proportional to the 
frequency uo, with dn/duo = NT. Any uncertainty An in 
the population measurement limits the attainable spec¬ 
troscopic sensitivity to AuoT > An/N. Correspondingly, 
approaching the Heisenberg limit requires single-particle- 
resolved state detection, which becomes increasingly dif¬ 
ficult at large atom number. Recent experiments have 
made progress in addressing this challenge mm, but 
not yet under the conditions required to generate highly 
entangled states. 

Theoretically, a quantum-enhanced measurement does 
not require directly detecting the entangled sensor state. 
Several proposals instead envision echo protocols [T4lfT6] 
in which a quantum system undergoes a unitary evolu¬ 
tion U into a non-classical state and, after subjecting 
this state to a perturbation, one attempts to reverse the 
evolution to the initial state by application of W. This 
approach in principle permits Heisenberg-limited mea¬ 
surements with an ensemble of spin-1/2 particles (or, 
equivalently, two-level atoms) by detecting the state of 
a single ancilla spin [[16] . 

In this Letter, we propose an echo protocol that en¬ 


ables spectroscopy near the Heisenberg limit with low- 
resolution state detection An ~ y/N. To generate en¬ 
tanglement, our method employs the global Ising inter¬ 
actions of the “one-axis twisting” Hamiltonian m , re¬ 
alizable with cold atoms mBnnEHE], trapped ions [23] 
EU and solid-state nuclear spins [26]. Switching the sign 
of the interaction after subjecting the system to a weak 
perturbation amplifies the perturbation into a larger spin 
rotation that is easily detected. We analyze the perfor¬ 
mance including dissipation in two atomic implementa¬ 
tions, employing interactions mediated either by light in 
an optical cavity T9j [21] or by Rydberg dressing [27]. In 
each case, the twisting echo enables precision far beyond 
the standard quantum limit with detection noise larger 
than the quantum noise of an unentangled state. 

The one-axis twisting Hamiltonian H tw [ s t = x$z de- 
scribes internal-state-dependent interactions in a collec¬ 
tion of N two-level atoms, which we represent in terms 
of spin-1/2 operators s$ by a collective spin S = J2iLi s i- 
The dynamical effect of iL tw ist is to generate a spin pre¬ 
cession about the z-axis at a rate proportional to S z (Fig. 
[l^i-b). For spins initially polarized along x, the lowest- 
order effect of i7 tw ist is squeezing PH- At longer interac¬ 
tion time, i^twist produces non-Gaussian states, including 
oversqueezed states and ultimately a maximally entan¬ 
gled GHZ state at = 7r/2 [24] . While the GHZ state 
enables Heisenberg-limited measurements in few-particle 
systems with highly coherent interactions, we will show 
that the twisting echo protocol attains a Heisenberg scal¬ 
ing AuoT oc 1/N at significantly shorter evolution time. 

The squeezed and oversqueezed states generated by 
#twist are highly sensitive to spin rotations fZ y ((f>) = 
e ~ i (f)Sy a pout the y axis. Indicative of this sensitivity 
are the increased quantum fluctuations A S y in Fig. [Tjo, 
which lower the quantum Cramer-Rao bound on the un¬ 
certainty A<f> > 1/(2A S y ) [28]. We therefore present a 
protocol for measuring fZ y (c £), bearing in mind that a 
compound sequence fZ y ((j)) = 1Z x {—tt/2)1Z z ((I))71 x (7t/2) 
then allows for measuring a precession </> = uoT about 
z, much as in Ramsey spectroscopy with squeezed states 
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The twisting echo protocol is shown in Fig. [lj where we 
assume unitary dynamics. An ensemble is initialized in 
the coherent spin state (CSS) |x) satisfying S x \x) = *S|x) 
(Fig. ^). Applying H twist (x) = xS 2 z for a time t yields 
the entangled state |^ e ) = U |x), where U = e~ lxS (Fig. 
[TJd) . To detect a rotation \x/j e ) 7Z y ((j))\i/j e ) by a small 
angle 0, we attempt to undo the twisting by applying 
-^twist(—x)• For 0 = 0, the final state U^lZ y {(j))U |x) is 
identical to the original CSS. However, a non-zero angle 
0 (Fig. |Tp biases the ^-dependent spin precession to 
produce a large final value of (S y ) (Fig. E*)- Measuring 
S y , by rotating the state and then detecting the popula¬ 
tion difference — n^, provides a sensitive estimate of 
0 . 
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FIG. 1. Twisting echo for entanglement-enhanced mea¬ 
surement. Top row: The initial CSS |x) (a) evolves under 
Htwist(x) an oversc l ueeze( l state |0 e ) (b). To detect a ro¬ 
tation of 10 e ) about y by a small angle 0 (b—>>c), we amplify 
the perturbation into a large displacement {S$) = GScf) by ap¬ 
plying H tw ist(—x) (c—kl). Illustrated are Wigner quasiproba¬ 
bility distributions for 2 S = 30 atoms, with 0 = 1/S. Bottom 
row: Cartoon depiction of the same steps, with blue flow lines 
indicating twisting and untwisting. 

The angular sensitivity is given by 

A<p=[AS$/d 4 ,(S$)]^ 0 , (1) 

where (Sfi) and A Sfi represent the mean and standard 
deviation of S y after the echo, and d$ = d/d(j). The 
standard deviation for no rotation is ideally that of the 
initial CSS, A S^ =0 = AS css = y/S/ 2. To evaluate the 
denominator of Eq. [l] we expand 

(S$) = {±\U ] e~ i(t)Sy US y U ] e i(t)Sy U\^) 

= ict>(x.\[S y ,UiS y U}\i)+0(<f> 2 ). (2) 

to lowest order in 0. We express S y = (S+ — S-)/(2i) in 
terms of raising and lowering operators S± and simplify 

U'S±U = e~ ixtS *S±e ixtS * = S ± e ixt ^ ±2Sz ^ (3) 

to evaluate Eq. [2] using the generating functions in Ref. 
[29] . We thus arrive at a dependence 

[Wv>]*=,o = S ( 2S ~ V 8 ™ (^) cos2S “ 2 (^) ( 4 ) 

of the final spin orientation on the perturbation 0, where 
we have introduced the ‘twisting strength’ Q = 2S'yt. 


The resulting metrological gain 1 /[N (A0) 2 ] is plotted in 
Fig. [2^i as a function of Q for N = 10 3 atoms. At the 
optimal twisting strength Q op t = 2 S arccot(\/2S' — 2) « 
y/~N for N 1, the echo protocol yields an angular sen¬ 
sitivity 

A0 min = y/e/N. (5) 

This sensitivity is very near the Heisenberg limit, despite 
a ^ v^-times shorter twisting evolution Q opt than re¬ 
quired to reach a GHZ state (Qghz in Fig. |2^i). The 
entangled state |0 e ) at Q 0 pt is oversqueezed (Fig. [lb), 
allowing the echo to surpass the sensitivity A0 oc 1 /TV 3 / 6 
attainable by spin squeezing m under H tw[st . 

The twisting echo is highly robust against detection 
noise (Fig. [2)3), as the “untwisting” amplifies the spin 
rotation signal by a factor of G = d(S^) / d(S 0) < y/N 
(Fig. -d). Concomitantly, the quantum noise returns 
to the CSS level, so that adding Gaussian detection 
noise AS meas = pAScss results in an angular sensitiv- 
ity A (j> = y/1 + p 2 A^ min . Thus, even a measurement 
that barely resolves a CSS, with atom number resolu¬ 
tion An = 2AS' meas ~ V~N, permits a sensitivity near 
the Heisenberg limit. By contrast, measurement noise 
significantly degrades the sensitivity attainable by direct 
detection of non-Gaussian states: already at single-atom 
resolution, the twisting echo outperforms direct detection 
of a GHZ state (Fig. [2)3). 

In practice, the sensitivity A0 may be degraded by 
imperfect coherence of the one-axis twisting evolution. 
To show that the twisting echo can provide a significant 
benefit in realistic metrological scenarios, we analyze the 
limitations due to dissipation in two implementations de¬ 
signed to enhance atomic clocks: the method of cavity 
feedback dynamics m demonstrated in Refs. H ED]; 
and the Rydberg dressing scheme proposed in Ref. m- 

Cavity-Mediated Interactions — The scheme for one- 
axis twisting by light-mediated interactions |T9ti2T] is 
shown in Fig. [3^i. Atoms in hyperfine states |t), |0) are 
coupled to an optical resonator mode with vacuum Rabi 
frequency 2p, at large detunings ±A from transitions to 
an excited state \e). The dispersive atom-light interac¬ 
tion shifts the cavity resonance frequency in proportion 
to S z , with duj c /dS z = <£>ft/2, where <£> = Ag 2 / An. Thus, 
driving the cavity at a detuning 5 C from the bare-cavity 
resonance results in an -dependent intracavity power. 
The latter acts back on the atomic levels via the a.c. 
Stark shift, yielding an -dependent spin precession. 
For small cavity shifts S c , the spin preces¬ 

sion rate depends linearly on S z , yielding one-axis twist¬ 
ing dynamics. The sign of the twisting is controlled by 
5 C , while the strength depends on the average number 
of photons p in the coherent field incident on the cavity: 
Q = 2 Sp& 2 d/ (1 + d 2 ) 2 , where d = 2 S c /n. 

The light-induced twisting is accompanied by fluctu¬ 
ations in the phase of the collective spin due to pho¬ 
ton shot noise. These fluctuations are described by a 
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FIG. 2. Metrological gain of the twisting echo with 
N = 10 3 atoms. Horizontal lines indicate the standard quan¬ 
tum limit (black) and Heisenberg limit (red), (a) Metrologi¬ 
cal gain vs. twisting strength Q for the unitary twisting echo 
(solid blue), compared to spin squeezing (dot-dashed green) 
and the quantum Cramer-Rao bound (QCRB) on phase sen¬ 
sitivity (dashed orange). The twisting echo nearly follows the 
QCRB to its plateau at Q ~ VN ; only at a much longer time 
Qghz = Ntt/2 does the QCRB increase by 3 dB to reach 
the Heisenberg limit, (b) Metrological gain vs. measurement 
uncertainty An for echo with twisting strength Q 0 pt (solid 
blue) or Q sq (dotted blue), compared to direct detection of 
the squeezed state at Q sq (dot-dashed green). Dashed orange 
curve shows Cramer-Rao bound for estimating cj) in the GHZ 
state lZ y ((/))( |y) + |— y)) using projective measurements with 
uncertainty An. 


Lindblad operator L = y/jS z in the master equation 
for the density matrix of the spin subsystem, where 
7 = 2 x/d [30]. Physically, yt represents the average 
number of photons rescattered into the cavity per atom 
while twisting. The leakage of these photons from the 
cavity in principle enables a measurement of S z , whose 
backaction is the dephasing described by a decay in 
((S± -0 ) 2 ) = e _47t ((S± _0 ’^ _0 ) 2 ) after twisting and un¬ 
twisting. Accordingly, to lowest order in yt <C 1, for 
N 1 , the variance of S y grows to 

(A S$= 0 ) 2 /AS'css 2 « 1 + 4 S^t = 1 + 4 Q/d. (6) 

Thus, cavity decay increases Ac/) by a factor yj 1 + 4 Q/d 
compared with the unitary case. 

The phase broadening can be made arbitarily small at 
large detuning S c at the price of increased decoherence 
from spontaneous emission. The latter occurs at a rate 
r sc = x(l/d + d)/(2rj) per atom, where 7 = 4g 2 /(ftT) is 
the single-atom cooperativity and T is the excited-state 
linewidth. Assuming each spontaneous emission event 
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FIG. 3. Experimental schemes for entanglement-enhanced 
measurement, (a) Optical cavity with drive laser at detuning 
S c from cavity resonance uj c and detunings ±A from atomic 
transitions |f) — > |e), \\) — >• |e). (b) Optical-lattice clock 
with metastable state |t) coupled to a Rydberg state | R) by a 
laser with Rabi frequency Q, producing the two-body poten¬ 
tial V (r) [27]. Both schemes yield one-axis twisting dynamics, 
with the sign of Htwist dictated by the sign of S c ( S R ) in the 
cavity-QED (Rydberg) system. 


has a probability r of flipping a spin, the average value of 
S z while “untwisting” differs from that during “twisting” 
by a root-mean-square value m A5f » y4r5r sc i/3 = 
\JrQ(l/d + d)/(‘3ri). Such a change has the same effect 
on the final signal as a perturbation A0 SC = A Sl c /S, and 
thus contributes to the uncertainty in measuring </>. 

To calculate the phase sensitivity of the dissipative 
echo, we first express the normalized phase variance 
ctq = 2 S (A^)q for the unitary case in terms of the 
twisting strength Q. From Eq. [4] using the approxi¬ 
mation cos 2S ~ 2 (xt ) ~ e ~ s (xt) f or yi we obtain 

ctq ~ e® /( 2S ^IQ 2 . The total phase variance including 
cavity decay and spontaneous emission is then given by 
a 2 = 25 (A (f>f = erg + cTj iss , where 

^diss = 4 e^/^/{Qd) + 2 rQ {dr 1 + d) /(3 Sr,) (7) 

represents the noise added by dissipation. This noise will 
reduce the optimal twisting strength below Q 0pt = \/2A, 
so that e® 2 /( 2S ^ ~ 1. The dissipative contribution cr| iss 
is then minimized by choosing Qd ~ yj6Srj/r and large 
detuning d^> 1, yielding a total variance 

a 2 « r( 1 + d 2 )/(6Sr]) + \/32r(l + l/d 2 )/(3Sr]). (8) 

At large collective cooperativity 2 Srj 1, the normalized 
phase variance a 2 depends only weakly on detuning (Fig. 
ka.i) for 1 d <C {/Srj/r , where we obtain a 2 ~ cr| iss ~ 
yj32r/(3Sr ]). We plot the metrological gain a~ 2 as a 
function of N and r] in Fig. [4^i.ii. 

The scaling of the metrological gain with collective 
cooperativity is just as in spin squeezing by quantum 
non-demolition measurement [33], but no coherence¬ 
preserving measurement is required. Moreover, the mea¬ 
surement only requires a resolution on the order of the 
width A > AS'css of the broadened final state. The 
benefits of the echo for state detection can be combined 
with schemes for reducing dissipation [13, HO [35] to 
achieve even higher sensitivity. 
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FIG. 4. Dissipative twisting echo, (a) Cavity-mediated 
twisting, with r = 1/2. (i) Metrological gain —101og 10 cr 2 

vs. laser-cavity detuning d = 2 8 c /k for N = 10 5 atoms and 
cooperativity rj = 0.1 (blue), 1 (orange), 10 (green), (ii) 
Metrological gain vs. N and 77 . (b) Rydberg-mediated twist¬ 
ing in Sr for n — 70, 2 a = 800 nm, e = 0.1. (i) Limitations 
on the detuning S R , which must lie below the blue line for 
all atoms to fit inside the interaction sphere, but above the 
dotted orange line to avoid spontaneous emission, (ii) Metro- 
logical gain vs. N in Rydberg scheme limited by spontaneous 
emission (orange), compared to idealized case without spon¬ 
taneous emission (dashed blue) and cavity-mediated twisting 
with 77 = 10 (dotted green). Orange band represents a range 
of Cq coefficients 10 10 < Cq < 10 11 calculated for principal 
quantum numbers 60 < n < 70 m- 


per atom, where T is the Rydberg state linewidth. To 
fit all atoms inside a sphere of diameter L in a three- 
dimensional lattice with spacing a, we must restrict the 
detuning to 5 R < TCq/(2 9 N 2 ), where Cq = Cq/(TcP). 
To conservatively estimate the metrological gain attain¬ 
able under this restriction, without modeling the effects 
of decay to other Rydberg states, we restrict the twisting¬ 
untwisting evolution to the first spontaneous emission 
event: 2NT sc t < 1. Reaching the optimum oversqueezed 
state at Qopt then requires a detuning S R > TN 3 ^ 2 /(2e). 
Both conditions on S R , plotted in Fig. [4}o.i, can be met 
simultaneously for up to N cr = (Cq6/2 3 ) 2 / 7 atoms. Here, 
Cq ~ 10 11 for a Rydberg state of principal quantum num¬ 
ber n ~ 70 [32: in a lattice of the “magic” wavelength 
2 a ~ 800 nm for the clock transition [36]. At a Rydberg- 
state population e = 0.1, the ideal phase sensitivity of 
Eq. [5] can then be reached with up to N cr ~ 150 atoms. 


We compare the predicted performance of the Rydberg 
and cavity schemes in Fig. .ii. For low atom number, 
the Rydberg dressing outperforms cavity-mediated inter¬ 
actions even at strong coupling 77 ^ 10. Yet whereas the 
cavity echo improves monotonically with TV, the Rydberg 
echo reaches an optimum at the critical atom number 
TV cr , above which the coherent evolution time must be 
reduced to extend the interaction range. Even with only 
TV cr ^150 atoms, the twisting echo matches the phase 
sensitivity of ^ 10 4 unentangled atoms. The method 
could thus benefit atomic clocks employing asynchronous 
interrogation of many small sub-ensembles [37, EH] • 


Rydberg Dressing —To approach the ideal unitary echo, 
we consider implementing H tw [ st by Rydberg dressing 
(Fig. [3 Jd) in a strontium optical-lattice clock [27] . Here, 
the pseudo-spin states are the singlet ground state |^) = 
I 1 So) and metastable triplet state |t) = | 3 Po)- A laser 
with Rabi frequency Q is detuned by S R from the |t) 

\R) transition, where \R) = \n 3 Si) represents a Ryd¬ 
berg state of principal quantum number n. The coupling 
to | R) induces a two-body potential V(r) that is nearly 
constant over a distance L = ^[Cq/( 2S r )] 1 ^ 6 (Fig. pp), 
enabling all-to-all interactions among TV atoms confined 
within a region of diameter L. In the limit of weak dress¬ 
ing, where the probability e = TVff2 2 /(4J 2 ) for even a 
single one of the TVf ~ TV/2 spin-up atoms to be excited 
to | R) is small (e 1), the result is a one-axis twisting 
Hamiltonian with interaction strength x = 12 4 /(16J 3 ) 
m- The sign of the interaction is controlled by the de¬ 
tuning , while the strength \x\ = 6 2 Cq/(2 5 N 2 L 6 ) is 
highest for small, dense ensembles. 

Obtaining maximally coherent all-to-all interactions 
requires careful choice of the detuning from the Ryd¬ 
berg state. The detuning S R controls both the interaction 
range L and the ratio x/^sc oc S R between the spin-spin 
coupling and the spontaneous emission rate T sc = eT/N 


We have presented a protocol that amplifies the effect 
of a phase rotation on an entangled state to enhance sig¬ 
nal readout. By transferring the phase information to 
the average displacement of a near-Gaussian state, the 
twisting echo attains a Heisenberg scaling in sensitiv¬ 
ity without single-particle resolution, and eliminates the 
need for Bayesian estimation methods for non-Gaussian 
states. Our approach can guide the design of new experi¬ 
ments by alleviating the need to simultaneously optimize 
coherence of interactions and fidelity of state detection. 
The protocol is adaptable to a wide range of systems, in¬ 
cluding ones where the sign of the interaction is fixed. For 
example, spin rotations can convert one-axis twisting to a 
two-axis counter-twisting Hamiltonian H ct oc S 2 — S 2 |55J 
and can switch the sign of H ct to exchange the squeezed 
and amplified quadratures. Future work might explore 
extensions to richer many-body systems featuring finite- 
range interactions or chaotic dynamics. 
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Approaching the Heisenberg Limit Without Single-Particle Detection: 

Supplemental Material 


In this supplement, we derive the effective one-axis twisting dynamics of an atomic ensemble coupled to a driven 
optical cavity [1-3]. This system has been analyzed elsewhere, but we know of no comprehensive treatment using the 
Lindblad master equation formalism. We also establish an equivalence between the twisting dynamics of the driven 
cavity scheme and the twisting dynamics of a four-photon Raman scheme proposed by Sprensen and Mplmer [4]. 


PHYSICAL SYSTEM 

The model system consists of N three-level atoms coupled to single cavity mode with annihilation operator a. The 
total Hamiltonian is: 


H — H c T H a T -Hint T H^ r [ ve 


(SI) 


where H c ,H a are the cavity and atomic Hamiltonians, respectively: 

H c := uo c a) a 

N 

H a := ^2 [2A |t) (t| i + (^c + A) \e) (e\ •]. (S2) 

i— 1 

The two atomic ground-state levels \l) i , |t)* are hyperfine sublevels separated by an energy difference 2A and the 
excited state \e) i is accessible from the ground states by optical transitions (see Fig. 3a in the main text). The cavity 
frequency uj c is detuned an equal amount ±A from the ground states |t)i > It)*- 

The dipole interaction term H int describes the coupling between the atoms and the cavity, 

N 

Hint ■= Y \-9 a \ e ) (It + ga \e) + h.c.], (S3) 

i= 1 


which excites the atom to \e) i by annihilating a cavity photon with Rabi frequency 2 g. Finally, the drive term Hd r i ve 
populates the cavity with photons from an external coherent driving field /3 at frequency uji: 


Hdrive .— 




a e 


—iuoit 


+ P a) 


jiuit 1 


(S4) 


where k is the cavity linewidth. In general, the drive laser is detuned from cavity resonance by an amount S c := uji — uo c . 

Including dissipative processes, the dynamics are described by a Lindblad master equation for the density matrix 
P- 


P = 


P] + 2 ( 


Up, L] 


L 




(S5) 


where Lindblad relaxation operators Lj describe the various decay processes. We introduce a relaxation operator L 
to describe the loss of photons from the cavity, 


L = yjh a, (S6) 

and 2 N relaxation operators for the spontaneous emission processes, 

Li, a = >/VW){e\ i , (S7) 

where Y is the spontaneous emission rate and a =t, i labels the atomic state after emission. 

We are interested in the effective dynamics of the ground states 1^ , |t)^. This is a large Hilbert space, but we can 
restrict the dynamics to a much smaller subspace by noting that the Hamiltonian couples the atoms to the cavity 
symmetrically. Therefore the atomic ensemble can simply be viewed as a single spin-A/2 that lives in the (N + 1)- 
dimensional symmetric Dicke subspace. Unfortunately, the spontaneous emission operators L* ?cr break this symmetry, 
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so in this analysis we will ignore them by assuming that the twisting dynamics occur before we lose a significant 
fraction of photons through emission. (The limitations imposed by spontaneous emission on the twisting evolution 
are analyzed heuristically in the main text.) 

The goal is to derive the effective dynamics of this spin-iV/2 in the presence of the driven cavity. We do this by 
adiabatically eliminating first the excited states \e) i and then the cavity a from the dynamics, leaving an effective 
Hamiltonian describing only the dynamics of the spin. To accomplish this, we follow the procedure of Reiter and 
Sorenson [5]. This procedure relies on two assumptions: 

1. The dynamics of the ground states , |t)i are slow compared to the dynamics of the rest of the system, and 

2. The ground states ||) •, \ J [) i are only perturbatively affected by the fast dynamics of the rest of the system. 

The first assumption allows us to solve for the fast dynamics in the presence of essentially stationary ground states, 
while the second ensures that the effective ground state dynamics remain slow even when coupled to the rest of the 
system. We refer the reader to the original paper for the full discussion. 


ADIABATIC ELIMINATION OF |e). 

We first eliminate the excited states \e) i using the Reiter-Sorenson procedure. To do so we must make two assump¬ 
tions: 

1. y/R\p\ ^yJa^a} <C A (i.e. excited state dynamics are fast compared to other timescales) 

2. g (ata) <C A (i.e. coupling to the excited states is perturbative) 

The second assumption is necessary in order to ignore the ground-state driving terms ~ e~ l6ct compared to 

the ground-state splitting 2A. 

Ignoring overall energy shifts, the result of this procedure is: 

H eS = K + §kS z / 2) a) a + 2A S z + [/ 3*a e iu)lt + pa) e“ iWit ] 

L e ff = \[k a (S8) 

where := 4g 2 /A (assuming T « A) and S z := \ JT a l z = \ JT [|^) (^|. — ||) (||.]. We see that the effect of the 
atomic ensemble is to shift the cavity resonance by QkS z / 2. In particular, we note that this cavity shift depends on 
the ^-projection S z of the ensemble spin. 

Passing to a rotating frame with respect to: 


Hq = 2A S z + uoia^a 


(S9) 


we obtain a time-independent effective Hamiltonian: 

Kft = (Sc + *kS z /2) ^ a + [p*a + pa)] 

L'efi = \/~k a (S10) 


MACROSCOPIC POPULATION OF CAVITY 

Due to the drive laser, the cavity is highly populated with photons, so we are not justified in simply adiabatically 
eliminating the occupied cavity mode a. Formally, the existence of a strong drive laser and populated cavity breaks 
the second assumption required for adiabatic elimination because the coupling terms ^ are not perturbative. 

However, we can always decompose the cavity field a into a large classical field a accompanied by small quantum 
fluctuations c: 


a = a + c. 


(Sll) 
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(Alternatively, this can be viewed as a displacement of the cavity operator a by an amount a.) We take the classical 
field a to be the steady-state population that would be induced by the driving field (3 in the absence of atoms: 

a = 13 \[^2 5 C + in/2 ^ S12 ^ 

Then the operator c describes the fluctuations in the cavity field due to its interaction with the atoms, which can be 
treated perturbatively. Plugging Eq. S12 into and L' eE we obtain: 

H" s = (S c + $kS z / 2) c f c + a*<f>nS z /2 c + a<f>K,S z /2 c f + <S>kS z /2 \a\ 2 

Lgg = \[k c (S13) 

where we have redefined the Hamiltionian to absorb the constant term that appears in L' eE . 


ADIABATIC ELIMINATION OF CAVITY FLUCTUATIONS 


We now eliminate the cavity fluctuations c. Applying the Reiter-Sorenson procedure once again, we obtain the 
result: 


-f^spin — 


^spm 


-S c + $kS z /2 


(~5 C + $kS x /2Y + k 2 /4 


$ 2 k 2 \a\ /4 S 2 


|a| /2 S z 


(S14) 


(~S C + §kS z /2Y + «2/4 

where we have made the following assumption: 

3. &kS z \a\ /2 <C \S C \ + $k,S z /2 (i.e. coupling to excited states is perturbative) 

The spin Hamiltonian simplifies further if in addition we assume that the shift &k,S z /2 in cavity frequency is small 

4. $«S*/2< \S C \ 

In this case we obtain: 


#sni 


spin 


xs 2 . 


-^spin — V^Td S z 


where we have defined: 


d := 2 8 C /k 

d 
X : 


(i + <py 


d P$ 2 


dt 


(S15) 


(S16) 


and dp/dt = \(3\ 2 is the rate at which coherent photons are incident on the cavity. During this analysis, we ignored 
the effect of spontaneous emission. We can, however, relate the rate of spontaneously emitted photons per atom to 
the interaction strength x as: 


( 9 _\ 

va; 


4- l«rr = 


Ixl 1 + d? 

2 r] d 


(S17) 


where rj := 4 g 2 /(kT) is the single-atom cooperativity. We use Eq. S17 to estimate the effects of spontaneous emission 
in the main text. (While this estimate does not treat the dynamical effects of spontaneous emission, it does give a 
sense of the timescale ~ 1/(7VT SC ) over which we can reasonably ignore these effects.) 

We reiterate the four assumptions that were made during this derivation, written in terms of the parameters 
d, 4>, dp/dt: 


i 

-L • 


dt V 1+d 2 


< A 
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3. \j dv ^ 1 j^p &S Z « d (using assumption 4 to ignore &S Z ) 

4. < d 


In our system, where the atoms are initialized in a coherent spin state in the x-y plane, assumption 4 is satisfied by 
keeping 4>V^ <C d, and assumptions 1-3 can always be satisfied by reducing the rate dp/dt of photons incident on 
the cavity. 

As a final note, it turns out that the effective spin dynamics in equation (S15) are equivalent to the dynamics 
one obtains from another one-axis twisting proposal that at first glance appears to be quite different. Sprenson and 
Mplmer proposed a four-photon Raman scheme to generate a one-axis twisting Hamiltonian by directly driving an 
atomic ensemble (e.g., using beams transverse to the cavity); the atoms communicate with one another by exchanging 
photons via a weakly occupied cavity mode [4]. The dynamics of this scheme are described by the effective operators: 

d^spin, SM = £ 

A S pin, SM = VW/d s x (S18) 

where we’ve defined 


£:= 


d 

1 +d 2 3 4 5 



(S19) 


and Qr is the Rabi frequency of the transverse drive lasers. In this scheme, the quantities A, S c are the single-photon 
detuning and two-photon detuning, respectively (other parameters, such as the cavity coupling g , are common to both 
schemes). The rate of spontaneous emission per atom in this scheme is: 


r sc , sm — 2 


(n R \\_ jgj ! + d 2 

V2Ay 2r) d 


(S20) 


Comparing equation (S15) to (S18) and equation (S17) to (S20) we see that, up to a rotation of spin axes S x —>> S z , 
the two schemes produce precisely the same dynamics. 

Although the two schemes are equivalent, the Sprenson-Mplmer scheme may offer an advantage for phase detection 
in that it twists about the S x axis. Thus, in an echo sequence that begins with a CSS along y, the twisting generates 
a state that is sensitive to a rotation about S z , which the untwisting amplifies into a large signal in (S z ) that can 
directly be detected in the populations of states |t), ||). By contrast, in the driven cavity scheme one must apply two 
extra rotations to perform an amplified phase measurement. 
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